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Abstract 

This chapter is an attempt to present a mathematical theory of compound fractional Poisson 
processes. The chapter begins with the characterization of a well-known Levy process: The com- 
pound Poisson process. The semi-Markov extension of the compound Poisson process naturally 
leads to the compound fractional Poisson process, where the Poisson counting process is replaced 
by the Mittag-Leffler counting process also known as fractional Poisson process. This process is no 
longer Markovian and Levy. However, several analytical results are available and some of them are 
discussed here. The functional limit of the compound Poisson process is an a-stable Levy process, 
whereas in the case of the compound fractional Poisson process, one gets an a-stable Levy process 
subordinated to the fractional Poisson process. 

PACS numbers: 
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I. INTRODUCTORY NOTES 



This chapter is an attempt to present a mathematical theory of compound fractional 
Poisson processes. It is not completely self-contained. The proofs of some statements can 
be found in widely available textbooks or papers. In several cases, freely downloadable 
versions of these papers can be easily retrieved. 

The chapter begins with the characterization of a well-known Levy process: The com- 
pound Poisson process. This process is extensively discussed in the classical books by Feller 
I| and de Finetti fl. 

The semi-Markov extension of the compound Poisson process naturally leads to the com- 
pound fractional Poisson process, where the Poisson counting process is replaced by the 
Mittag-Leffler counting process also called fractional Poisson process (3), 0]. This process is 
no longer Markovian and Levy. However, several analytical results are available and some 
of them are discussed below. 

The functional limit of the compound Poisson process is an a-stable Levy process, whereas 
in the case of the compound fractional Poisson process, one gets an a-stable Levy process 
subordinated to the fractional Poisson process. 

I became interested in these processes as possible models for tick-by-tick financial data. 
The main results obtained by my co-workers and myself are described in a review paper for 
physicists Q|. 



The reader interested in Monte Carlo simulations can consult two recent papers 
where algorithms are presented to simulate the fractional compound Poisson process. 



II. COMPOUND POISSON PROCESS AND GENERALIZATIONS 

Let {X i \ c *L 1 be a sequence of independent and identically distributed (i.i.d.) real-valued 
random variables with cumulative distribution function Fx{x), and let N(t), t > denote 
the Poisson process. Further assume that the i.i.d. sequence and the Poisson process are 
independent. We have the following 

Definition II. 1 (Compound Poisson process). The stochastic process 

N(t) 

Y(t) = Y,Xt (!) 
i=i 
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is called compound Poisson process. 

Here, we shall consider the one- dimensional case only. The extension of many results to 
M. d is often straightforward. The compound Poisson process can be seen as a random walk 
subordinated to a Poisson process; in other words, it is a random sum of independent and 
identically distributed random variables. It turns out that the compound Poisson process is 
a Levy process. 

Definition II. 2 (Levy process). A stochastic process Y(t), t > with Y(0) = is a Levy 
process if the following three conditions are fulfilled 

1. (independent increments) if t\ < t 2 < ■■■ < t n , the increments Y(t 2 ) — 
Y(ti), . . . , Y(t n ) — Y(t n ^i) are independent random variables; 

2. (time-homogeneous increments) the law of the increment Y(t + At) — Y(t) does not 
depend on t; 

3. (stochastic continuity) Va > 0, one has that liniAt->o P(|^(t + At) — Y(t)\ > a) = 0. 

Loosely speaking, one can say that Levy processes are stochastic processes with stationary 
and independent increments. Due to Kolmogorov's extension theorem 8|, a stochastic pro- 
cess is characterized by its finite dimensional distributions. In the case of a Levy process, the 
knowledge of the law of Y(At) is sufficient to compute any finite dimensional distribution. 
Let us denote by fY(At){y, At) the probability density function of Y(At) 

f Ym (y, At) dy d ^ f F(Y(At) e dy). (2) 

As an example, suppose you want to know the joint density function /y(t 1 ),y(t 2 )(?/i, h; y 2 , t 2 ) 
defined as 



fY(t 1 ),Y(t 2 )(yi,ti;y 2 ,t 2 )dy 1 dy 2 = P(T(ti) e dy 1 ,Y(t 2 ) e dy 2 ). (3) 



This is given by 



fY(ti),Y(t2)(yh t i'iV2,h)dyidy 2 = P(Y(*i) G dy h Y(t 2 ) e dy 2 ) 



= F{Y(t 2 ) e dy 2 \Y{t x ) e dy 1 )¥(Y(t 1 ) e dy x ) 
F(Y(t 2 ) - F(tx) G d(y 2 - y x ))P(y(ti) G d Vl ) 

= fY(t 2 )-Y(t 1 )(y2 - yi,t 2 - t 1 )f Y (t 1 )(yi,ti)dyidy 2 , (4) 
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and this procedure can be used for any finite dimensional distribution. The extension the- 
orem shows the existence of a stochastic process given a suitable set of finite dimensional 
distributions obeying Komogorov's consistency conditions but not the uniqueness. 

Definition II. 3 (cadlag process). A stochastic process Y(t), t > is cadlag (continu a 
droite et limite a gauche) if its realizations are right- continuous with left limits. 

A cadlag stochastic process has realizations with jumps. Let i denote the epoch of a 
jump. Then, in a cadlag process, one has Y(t) = Y(t + ) = f lim i _ > j+ Y(t). 

Definition II. 4 (Modification of a process). A modification Z(t), t > 0, of a stochastic 
process Y{t), t > 0, is a stochastic process on the same probability space such that ¥(Z(t) = 
Y(t)) = 1. 

Theorem II. 1. Every Levy process has a unique cadlag modification. 

Proof. For a proof of this theorem one can see the first chapter of the book by Sato □ 

The following theorem gives a nice characterization of compound Poisson processes. 

Theorem II. 2. Y(t) is a compound Poisson process if and only if it is a Levy process and 
its realizations are piecewise constant cadlag functions. 

Proof. An accessible proof of this theorem can be found in the book by Cont and Tankov 
H. □ 

As a consequence of the above results, the compound Poisson process enjoys all the 
properties of Levy processes, including the Markov property. To show that a Levy process 
has the Markov property, some further definitions are necessary. 

Definition II. 5 (Filtration). A family Tt, t > of a -algebras is a filtration if it is non- 
decreasing, meaning that F s C JF t for < s < t. 

Definition II. 6 (Adapted process). A process Y(t), t > is said adapted to the filtration 
Tt, t > if it is Tt-measurable for t > 0. 

Definition II. 7 (Markov process with respect to a filtration). A process Y(t) is a Markov 
process with respect to the filtration T t , t > if it is adapted to T t , t > and (A C M.) 

F(Y(t) G A\F a ) = F(Y(t) G A\Y(s)). (5) 



Definition II. 8 (Natural filtration). The natural filtration for a stochastic process is the 
family of non- decreasing a-algebras generated by the process itself {a(X(s)), s e [0,t]} 7 
t > 0. 



Definition II. 9 (Markov process with respect to itself). A process Y(t), t > is a Markov 
process with respect to itself (or simply a Markov process) if it is a Markov process with 
respect to its natural filtration. 

The natural filtration {a(X(s)), s G [0, t]}, t > 0, is a formal way to characterize the 
history of the process up to time t. For a Markov process, the future values do not depend 
on the whole history, but only on the present value of the process. 

Definition 11.10 (Transition probability). Given the Markov process Y(t), t > 0, its tran- 
sition probability P(y, A, At, t) is defined as 



where id. 

Definition 11.11 (Homogeneous Markov process). A Markov process Y(t), t > is said 
(time) -homogeneous if its transition probability P(y, A, At,t) does not depend on t. 

Theorem II. 3. A Levy process is a time-homogeneous Markov process with transition prob- 
ability 



Proof. The Markov property is a consequence of the independence of increments. The 
following chain of equalities holds true 



P(y, A, At, t) = F(Y(t + At) e A\Y(t) = y), 



(6) 




(7) 



F(Y(t + At) e A\T t ) = F(Y(t + At) - Y(t) e A - Y(t)\F t ) 



F(Y(t + At) - Y(t) e A — Y(t)\Y(t)). 



(8) 



We further have 



F(Y(t + At)eA\Y(t)=y) = F(Y(At)eA-y)= [ f Y{At) (x - y, At) dx, (9) 



JxeA 



as a consequence of time homogeneity. 



□ 
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This result fully justifies the passages leading to equation fl4|. 

It turns out that fY(At){y, At) can be explicitly written for a compound Poisson process. 
Let F x (x) be the law of the jumps {Xi} and let A denote the parameter of the Poisson 
process, then we have the following 

Theorem II. 4. The cumulative distribution function of a compound Poisson process is given 
by 

(\ty 



F Y{t) (y,t) = e-*J2-^r F Y:(y)> ( 10 ) 

n=0 

where F Y ™(y) is the n-fold convolution of Fx(x). 

Proof. Starting from Y(0) = 0, at time t, there have been N(t) jumps, with N(t) assuming 
integer values starting from (N(t) = means no jumps up to time t). To fix the ideas, 
suppose that N(t) = n. Therefore, one has 

N(t) n 

Y{t) = Y,X l = Y, X ^ = Y n (11) 

i=l i=l 

and, in this case, 

F Yn {y) = <y)=F^2x i <y y j= F%(y). (12) 

For the Poisson process, one has 

P(n, t) = F(N(t) =n)= e~ xt ^^. (13) 

n\ 

Given the independence between N(t) and the XjS, one has that 

< V,N(t) =n) = P(n,t)F Y :(y) = e^^F*"^). (14) 
The events {Y n < y, N(t) = n} are mutually exclusive and exhaustive, meaning that 

{Y(t) <y} = U~ {Y n < y, N(t) = n}, (15) 

and that, for any m ^ n 

{Y m < V, N(t) =m}n {Y n < y, N(t) =n} = 0. (16) 
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Calculating the probability of the two sides in equation ( Fl5|) and using equation (TjrjJ and 
the axiom of infinite additivity 



F Y{t) (y,t) = F(Y(t) < y) = P(U~ {Y n < y,N(t) = n}) 

oo oo 

- xt E K -^r F r:(y) ( 17 ) 



n=0 n=0 

" (At)' 



n=0 

leads to the thesis. □ 

Remark II. 1 (Generality of Theorem III.4I) . Note that the above theorem is valid for any 
counting process N(t) in the following form 

Theorem II. 5. Let {X} c *l 1 be a sequence of i.i.d. real-valued random variables with cumu- 
lative distribution function Fx(x) and let N(t), t > denote a counting process independent 
of the previous sequence and such that the number of events in the interval [0, t] is a finite 
but arbitrary integer n = 0, 1, . . .. Let Y{t) denote the process 

N{t) 

Y{t) = Y J X l . (18) 

i=l 

Then if P(n,t) =^P(A r (t) = n), the cumulative distribution function ofY{t) is given by 

oo 

F Y{t) (y,t) = Y,P(n,t)F™(y), (19) 

n=0 

where F Y ™(y) is the n-fold convolution of Fx(x). 

Proof. The proof of this theorem runs exactly as the proof of theorem III.4I without specifying 
P(n,t). □ 

Theorem III. 51 will be useful in the next section when the Poisson process will be replaced 
by the fractional Poisson process. 

Remark II. 2 (The n = term). For n = 0, one assumes F Y ®{y) = 0(y) where 6(y) is the 
Heaviside function. Note that P(0, t) is nothing else but the survival function at y = of 
the counting process. Therefore, equation ( Tl9|) can be equivalently written as 

oo 

F Y{t) (y, t) = P{n, 0) e{y) + ^ P(n, t)F Y n n {y), (20) 

n=l 
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Remark II. 3 (Uniform convergence). The series (fTUj) and (|T9|) are uniformly convergent for 
y 7^ and for any value of t e (0, oo) (this statement can be proved applying Weierstrass 
M test). For y = there is a jump in the cumulative distribution function of amplitude 
P(0,t). 

Example II. 1 (The normal compound Poisson process). As an example of compound Pois- 
son process, consider the case in which Xi ~ A/"(/i,a 2 ), so that the cumulative distribution 
function is 

F x (x) = $(x|/i,a 2 ) = ~ (l + erf (^=)) , (21) 

where 

erf (x) = -= [ e~" 2 du (22) 



o 



7T „ 

2\ 



is the error function. In this case, the convolution F Y ™(y) is given by Q(y\rifi, na ) and one 
finds 

F Y{t) (y,t) =e- xt f^ { -^f$(y\n^na 2 ). (23) 

n=0 

Corollary II. 6. In the same hypotheses as in Theorem UI.b^ the probability density fv(t)(y, t) 
of the process Y(t) is given by 

oo 

f Y{t) (y, t) = P(0, t) 8{y) + p in, t)f£(y), (24) 

n=l 

where fy™(y) is the n-fold convolution of the probability density function fy n (y) = 
dFyM/dy- 

Proof. One has that /y(t)(y,t) = dF Y (t)(y,t)/dy; moreover, equation ([21]) is the formal 
term-by-term derivative of equation (1191) . If y ^ 0, there is no singular term and the series 
converges uniformly (/y™(y) is bounded and Weierstrass M test applies), therefore, for any 
y it converges to the derivative of FY(t){y,t). This is true also for y = for n > 1 and the 
jump in y = gives the singular term of weight P(0,t) (see equation (|20|) ). □ 

Remark II. 4 (Historical news and applications). The distribution in equation f|T9|) is also 



known as generalized Poisson law. : 
a famous work published in 1943 



his class of distributions was studied by W. Feller in 
111 ] . It is useful to quote an excerpt from Feller's paper, 



with notation adapted to the present paper. 
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The most frequently encountered application of the generalized Poisson distribu- 
tion is to problems of the following type. Consider independent random events 
for which the simple Poisson distribution may be assumed, such as: telephone 
calls, the occurrence of claims in an insurance company, fire accidents, sickness, 
and the like. With each event there may be associated a random variable X. 
Thus, in the above examples, X may represent the length of the ensuing conver- 
sation, the sum under risk, the damage, the cost (or length) of hospitalization, 
respectively. To mention an interesting example of a different type, A. Einstein 
Jr. [12] and G. Polya jl^, [l^ have studied a problem arising out of engineering 
practice connected with building of dams, where the events consists of the mo- 
tions of a stone at the bottom of a river; the variable X is the distance through 
which the stone moves down the river. 

Now, if F(x) is the cumulative distribution function of the variable X associated 
with a single event, then F* n (x) is the cumulative distribution function of the 
accumulated variable associated with n events. Hence the following equation 



oo r 

a r 



* nl x) 



nl 

n=0 



is the probability law of the sum of the variables (sum of the conversation times, 
total sum paid by the company, total damage, total distance travelled by the 
stone, etc.). 

In view of the above examples, it is not surprising that the law, or special cases 
of it, have been discovered, by various means and sometimes under disguised 
forms, by many authors. 

Indeed, the rediscovery and/or reinterpretation of equation ffl9|) went on also after Feller's 
paper. In physics, X is interpreted as the position of a walker on a lattice and N(t) is 
the number of walker jumps occured up to time t 15Ml9|. In finance, X is the tick-by-tick 
log-return for a stock and N(t) is the number of transactions up to time t j^. 

The application of Fourier and Laplace transforms to equation (12%] ) leads to an equation 
which is known as Montroll- Weiss equation in the physics literature [15|. For reasons which 
will become clear in the following, it can also be called semi-Markov renewal equation. Let 

r+oa 

/(«)= / e iKX f{x)dx 



5C 



denote the Fourier transform of f(x) and 

POO 

g(s)= / e- st g(t)dt 
Jo 

denote the Laplace transform of g(t), then the following theorem holds true. 

Theorem II. 7. (Montroll-Weiss equation) Let J denote the sojourn time of the Poisson 
process in N(t) = 0, with Fj(t) = 1 - e~ xt , fj(t) = \e~ xt and P(0,t) = 1 - Fj(t) = e~ xt . 
We have that: 

1. The Fourier-Laplace transform of the probability density f Y (t)(y,t) is given by 

f Y(t) ( K ,s) = ? (25) 

2. The probability density function f'Y(t){y,t) obeys the following integral equation 

f Y{t) (y,t) = P(0,t)5(y)+ dt'fj(t-t') dy'f x {y -y')f Y(t) {y'X). (26) 

JO J -oo 

In order to prove Theorem (111.71) . we need the following lemma. 

Lemma II. 8. Let T 1; T 2 , . . . , T n , . . . denote the epoch of the first, the second, . . ., the n-th, 
. . . event of a Poisson process, respectively. Let J — J\ — 7\ denote the initial sojourn time 
and, in general, let Ji = T{ — Tj_i be the i-th sojourn time. Then {Ji}^i is a sequence of 
i.i.d. random variables. 

Proof. The proof of this lemma can be derived combining lemma 2.1 and proposition 2.12 



in the book by Cont and Tankov 10]. □ 



It is now possible to prove the theorem 

Proof. (Theorem III. 711 Let us start from equation (124j) and compute its Fourier-Laplace 
transform. It is given by 

f Y(t) (K, s) = P(0, 3) + P(n, s) [&(«)] n . (27) 

71=1 

Now, we have that 

n 

T n = J2^ (28) 
i=i 
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is a sum of i.i.d. positive random variables and P(n,t) = f ¥(N(t) = n) meaning that 
there are n jumps up to t = t n and no jumps in t — t n . Therefore, from pure probabilistic 
considerations, one has that 

P(n,t) = P(Q,t-t n )*f Tn {t n ) (29) 
and, as a consequence of equation ( 128|) . one further has that 

fT n (Q = fT(tn). (30) 

Therefore, one can conclude that 

P(n,s) = P(0,s)[fj(s)\ . (31) 

Inserting this result in equation (|27|) . noting that P(0, s) = (1 — fj(s))/s, and summing the 
geometric series (fx and fj are probability densities) leads to equation ( |25j) : 



f Y(t) (K, S ) = P(0, S )+P(0,s)f^[fj( S )f X (K) 

n=l 

OD 

= P(0, S )J2[fj(s)fx(K 

n=0 

1 - f.i(s) 1 



S 1 " fj(s)fx(K) 

Equation (|25|) can be re-written as 



(32) 



/y(t)(«»«) = - — +/j(s)/x(«)/y (t )(/«,s); (33) 

Fourier-Laplace inverting and recalling the behaviour of convolutions under Fourier-Laplace 
transform, immediately leads to equation ( 12T)|) . □ 

Remark II. 5. Theorem ( 111.71) was proved in the hypothesis that N(t) is a Poisson process. 
In this case, one has P(0,s) = 1/(A + s) and fj(s) = A/(A + s) and equation ( 1251) becomes 

?y (0 («, S ) = — 4— — • (34) 

A - Xfx{K>) + s 

The inversion of the Laplace transform yields the characteristic function of the compound 
Poisson process 

f Y{t) (K,t) = E (e tKY{£) ) = e-W-fxW*. (35) 
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Remark II. 6. The proof of Theorem (1H.7|) does not depend on the specific form of P(0, t) 



and fj(t), provided that the positive random variables {J} C *L 1 are i.i.d.. Therefore, equations 
(|2^|) and (|2T)j) are true also in the case of general compound renewal processes starting from 
Y(0) = at time 0. 

III. COMPOUND FRACTIONAL POISSON PROCESSES 

Definition III.l (Renewal process). Let {7}°^ be a sequence of i.i.d. positive random 
variables interpreted as sojourn times between subsequent events arriving at random time. 
They define a renewal process whose epochs of renewal (time instants at which the events 
take place) are the random times {T}^ =0 defined by 

To = 0, 

n 

T n = ^2 Ji- ( 36 ) 

i=l 

The name renewal process is due to the fact that at any epoch of renewal, the process 
starts again from the beginning. 

Definition III. 2 (Counting process). Associated to any renewal process, there is a counting 
process N(t) defined as 

N(t) = max{n : T n < t} (37) 
that counts the number of events up to time t. 

Remark III.l. As mentioned in the previous section N(t) is the Poisson process if and 
only if J ~ exp(A). Incidentally, this is the only case of Levy and Markov counting process 



related to a renewal process (see Cinlar's book [20] for a proof of this statement). 



Remark III. 2. In this paper, we shall assume that the counting process has cadlag sample 
paths. This means that the realizations are represented by step functions. If t%. is the epoch 
of the k-th jump, we have N(t^) = k — 1 and N(t^) = k. 

Let {Xi}^ be a sequence of independent and identically distributed (i.i.d.) real-valued 
random variables and let N(t), t > denote the counting process. Further assume that the 
i.i.d. sequence and the counting process are independent. We have the following 
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Definition III. 3 (Compound renewal process). The stochastic process 



N(t) 

Y(t) = J2Xi (38) 
i=i 



is called compound renewal process. 



Remark III. 3. Again, here, it is assumed that the sample paths are represented by cadlag 
step functions. Compound renewal processes generalize compound Poisson processes and 
they are called continuous-time random walks in the physical literature. 

Remark III. 4. As compound renewal processes are just Markov chains (actually, random 
walks) subordinated to a counting process, their existence can be proved as a consequence 
of the existence of the corresponding discrete-time random walks and counting processes. 

In general, compound renewalpro cesses are non-Markovian, but they belong to the wider 



class of semi-Markov processes [7|, I20H23 ] . 



il_pro( 



Definition III. 4 (Markov renewal process). A Markov renewal process is a two-component 
Markov chain {Y n ,T n }^L , where Y n , n > is a Markov chain and T n , n > is the n-th 
epoch of a renewal process, homogeneous with respect to the second component and with 
transition probability defined by 

P(F n+1 G A, J n+1 < t\Y , . . . Y n , Ji, . . . , J n ) = F{Y n+1 G A, J n+1 < t\Y n ), (39) 

where A C R is a Borel set and, as usual, J n+ i = T n+1 — T n . 

Remark III. 5. In this paper, homogeneity with respect to the first component will be 
assumed as well. Namely, if Y n = x, the probability on the right-hand side of equation (1391) 
does not explicitly depend on n. 

Remark III. 6 (Semi-Markov kernel). The positive function Q(x, A,t) = P(y n+1 = y G 
A, J n +i < t\Y n = x), with x G R, A C R a Borel set, and t > is called semi-Markov kernel. 

Definition III. 5 (Semi-Markov process). Let N(t) denote the counting process defined as 
in equation l[3l\ ), the stochastic process Y(t) defined as 



Y(t) = Y m (40) 
is the semi-Markov process associated to the Markov renewal process Y n ,T n , n > 0. 
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Remark III. 7. In equation (l3"7|) . max is used instead of the more general sup as only 
processes with finite (but arbitrary) number of jumps in (0, t] are considered here. 

Theorem III.l. Compound renewal processes are semi-Markov processes with semi-Markov 
kernel given by 

Q(x,A,t)=P(x,A)Fj(t), (41) 

where P(x,A) is the Markov kernel (a.k.a. Markov transition function or transition proba- 
bility kernel) of the random walk 

P(x,A) = f nY n+1 G A\Y n = x), (42) 

and Fj(t) is the probability distribution function of sojourn times. Moreover, let fx(x) 
denote the probability density function of jumps, one has 

P(x,A) = [ fj(u)du, (43) 

JA-x 

where A — x is the set of values in A translated of x towards left. 

Proof. The compound renewal process is a semi-Markov process by construction, where the 
couple Y n ,T n , n > defining the corresponding Markov renewal process is made up of a 
random walk Y n , n > with Y = and a renewal process with epochs given by T n , n > 
with To = 0. Equation (HIT) is an immediate consequence of the independence between the 
random walk and the renewal process. Finally, equation fT4"31 is the standard Markov kernel 
of a random walk whose jumps are i.i.d. random variables with probability density function 
fx(x). □ 

Remark III. 8. As a direct consequence of the previous theorem, if the law of the couple 
X n , J n has a joint probability density function fx t j(x,t) = fx(x)fj(t), then one has 

P(F n+1 G A, J n+1 < t\Y n ) = Q{x,A,t) = P(x,A)Fj(t) = 

= [ f x (u)du [ fj(v)dv. (44) 

J A JO 

Theorem III. 2. (Semi-Markov renewal equation) The probability density function fY(t)(y, t) 
of a compound renewal process obeys the semi-Markov renewal equation 

pt r+oo 

f Y{t) (y,t)=P(0,t)8(y)+ dt'fj(t-f) dy'f x ( y - y ')f Y{t) ( y ',t'). (45) 

J0 J~oo 
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Proof. By definition, one has that 

P(Y(t) G dy\Y(0) = 0) = f Y (t){y,t)dy, (46) 

and that 

P(Y(t) G dy|Y(t') = y') = P(Y(t - f) G dy|Y(0) = </) = 
F(Y(t-t')-y' Edy\Y(0) = 0) = f Y{t) (y - y', t - t')dy, (47) 

because the increments in time and space are i.i.d. and hence homogeneous. From equation 
(14"4"|) . one further has 

P(Yi G dy, h G d*|Y = 0) = f x (y)fj(t)dydt. (48) 

Now, the probability in equation (l4"6"j) can be decomposed into two mutually exclusive parts, 
depending on the behaviour of the first interval 

P(Y(t) G dy\Y(0) = 0) = 
P(Y(t) G dy, Ji > t|y(0) = 0) + P(Y(t) G dy, Ji < t\Y(0) = 0). (49) 

The part with no jumps up to time t immediately gives 

P(Y(t) G dy, J x > t\Y(0) = 0) = P(0, t)S(y)dy, (50) 

whereas the part with jumps becomes 

F(Y(t) G dy, J x < t\Y(0) = 0) = 

/+oo />t 
/ P(Y(t) G ^lY(t') = y')P(Yi G Ji G dt') = 

-oo 

/ ~y',t- t')dyf x {y')fj{t')dy'dt' = 

x) JO 

/ + CXD /"t 
/ fY(t)(y-y\t-t')f x ( y ')fj(t')dy'dt' dy (51) 
_ -oo JO 

as a consequence of Bayes' formula and of equations ( 14T|) and f|48|) . A replacement of 
equations (1461) . ([50]) . (ISTjl into equation ([491 and a rearrangement of the convolution variables 
straightforwardly lead to the thesis ( 1451) . □ 

Remark III. 9. Note that the semi-Markov renewal equation ([45]) does coincide with the 
Montroll- Weiss equation f[26[) as anticipated. 



15 



Definition III. 6 (Mittag-Leffler renewal process). The sequence {J«}^ 1 of positive inde- 
pendent and identically distributed random variables whose complementary cumulative dis- 
tribution function Pp(Q, t) is given by 

Pp(0,t) = Epi-tP) (52) 

defines the so-called Mittag-Leffler renewal process. 

Remark III. 10. The one-parameter Mittag-Leffler function in fl52|) is a straightforward 
generalization of the exponential function. It is given by the following series 

oo n 

W = Sf(S+i)' (53) 

n=0 1 

where T(z) is Euler's Gamma function. The Mittag-Leffler function coincides with the 
exponential function for /3 = 1. The function Ep{—t^) is completely monotonic and it is 1 
for t = 0. Therefore, it is a legitimate survival function. 

Remark III. 11. The function Ep{—t^) is approximated by a stretched exponential for 
t -»■ 0: 

Ep{-tP) ~ 1 - T * ~ e -^ +1 ), for < t « 1, (54) 

and by a power-law for t — > oo: 

_, , sin(/37r) r(/3) . _ ,__ s 

E p {-tP)~ Kl ' y ;> fort»l. (55) 

Remark III. 12. For applications, it is often convenient to include a scale factor in the 
definition (152|) . so that one can write 

P(0,t) = E P (-(i/v)") . (56) 

As the scale factor can be introduced in different ways, the reader is warned to pay attention 
to its definition. The assumption j t = 1 made in (1521) is equivalent to a change of time unit. 

Theorem III. 3. (Mittag-Leffler counting process - fractional Poisson process) The counting 
process Np{t) associated to the renewal process defined by equation ( 152"]) has the following 
distribution 

Ppfat) = ¥(Np(t) = n) = z —E$\-iP), (57) 
where Ep(—t^) denotes the n-th derivative of E@(z) evaluated at the point z = —t 13 . 



16 



Proof. The Laplace transform of P(0,t) is given by 24 1 

as a consequence, the Laplace transform of the probability density function fj${t) = 
—dPp(Q,t)/dt is given by 

Ms) = (59) 
recalling equation (j3ip . one immediately has 

1 s?- 1 

Pa{n, s) = - s . (60) 

pv y (1 + s^) n 1 + v 7 

I — I 

Using equation (1.80) in Podlubny's book [24] for the inversion of the Laplace transform in 
f l60|) . one gets the thesis ([57]). □ 



Remark III. 13. The previous theorem was proved by Scalas et al. 25|, [261]. Notice that 
Ni(t) is the Poisson process with parameter A = 1. Recently, Meerschaert et al. |27J proved 
that the fractional Poisson process Np(t) coincides with the process defined by Ni(E(t)) 
where E(t) is the functional inverse of the standard /3-stable subordinator. The latter process 
was also known as fractal time Poisson process. This result unifies different approaches to 
fractional calculus 



28 



29J. 



Remark III. 14. For < (5 < 1, the fractional Poisson process is semi-Markov, but not 
Markovian and is not Levy. The process Np(t) is not Markovian as the only Markovian 
counting process is the Poisson process |20] . It is not Levy as its distribution is not infinitely 
divisible. 

Definition III. 7 (Compound fractional Poisson process). With the usual hypotheses, the 
process 

Nf>(t) 

Y&) = Y m) =Y. X * ( 61 ) 

i=i 

is called compound fractional Poisson process. 

Remark III. 15. The process Y\(t) coincides with the compound Poisson process of param- 
eter A = 1. 

Theorem III. 4. Let Yp(t) be a compound fractional Poisson process, then 
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1. its cumulative distribution function Fy (t)(y,t) is given by 



F m) (y,t) = E P {-i?)o{y) + ^r4 n) (-^)^:(i/) 



n=l 



(62) 



2. its probability density fy M(y,t) function is given by 



t ftr. 



n=l 



(63) 



3. its characteristic function /y»(t)(«, t) is given by 

f Y0it) (K,t)=E^y(f X (K)-i) 



(64) 



Proof. The first two equations (162]) and (163]) are a straightforward consequence of Theo- 
rem 111.51 Corollary III. 61 and Theorem IIII.3I Equation ( |64]) is the straightforward Fourier 
transform of ([6"3"]) . □ 

Remark III. 16. For < /3 < 1, the compound fractional Poisson process is not Markovian 
and not Levy (see Remark IIH. 141) . 



IV. LIMIT THEOREMS 



Definition IV. 1 (Space-time fractional diffusion equation). Let d a /d\x\ a denote the spatial 
non-local pseudo-differential operator whose Fourier transform is given by 

'd a f(x) 



T 



d\x\ 



-;« 



-N°7(«) 



(65) 



for x G (— oo, +oo) ; < a. < 2 and let /dt@ denote the time non-local pseudo-differential 
operator whose Laplace transform is given by 

'd?g{t)_ 



dt(> 



(66) 



for t > 0, < (3 < 1. Then the pseudo-differential equation 

d a u(x ) t) d^utx, t) 
d\x\ a = dtP 

is called space-time fractional differential equation. 



(67) 
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Remark IV. 1. The operator d a /d\x\ a is called Riesz derivative and is discussed by Saichev 



and Zaslavsky 



30L The operator d 13 /dt" is called Caputo derivative and was introduced by 



Caputo in 1967 31] as a regularization of the so-called Riemann-Liouville derivative. 



Theorem IV. 1. (Cauchy problem for the space-time fractional diffusion equation) Consider 
the following Cauchy problem for the space-time fractional diffusion equation (1671) 

d a u at p(x,t) d p u a ^(x,t) 



then the function 



where 



<9|x| a 



u a;/3 (x,t) 



me 



8{x) 



1 jxr 1 •'' 



t ft/a 



idle 



1 /- + 00 

W a , fi (u) = — / dKe- iKU Ep(-\K\ a ), 



(68) 



(69) 



(70) 



solves the Cauchy problem 



Proof. Taking into account the initial condition ( 1681) . as a consequence of the operator 
definition, for non-vanishing k and s, the Fourier-laplace transform of equation f )67|) is given 
by 

(71) 



leading to 



u a ,p{K, s) 



\n\ a + s^' 

Equation (1581) can be invoked for the inversion of the Laplace transform yielding 



u 



a ^(K,t)=Ep(-tP\K\ a ). 



Eventually, the inversion of the Fourier transform leads to the thesis. 



(72) 



(73) 



□ 



Remark IV. 2. The function defined by equations ( l69l) and ( |70l) is a probability density 
function. For — 1 and a = 2, it coincides with the Green function for the ordinary 
diffusion equation. The case /3 = 1 and < a < 2 gives the Green function and the 
;ransition probability density for the symmetric and isotropic a-stable Levy process L a (t) 



33]. 
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Theorem IV. 2. Let Y a ^(t) be a compound fractional Poisson process and let h and r be 
two scaling factors such that 

Y n (h) = hX 1 + . . . + hX n (74) 
T n (r) = rJi + . . . + rJ n , (75) 

and 

lim % = 1, (76) 

TOi/i < a < 2 and < /3 < 1. To clarify the role of the parameter a, further assume that, 
for h — > 0, one has 

f x (h K )~l-h a \K\ a , (77) 

then, for h, r — > with h a /r^ — » 1, f hYafj ( r t)(x,t) weakly converges to u a) p(x, t), the Green 
function of the fractional diffusion equation. 

Proof. In order to prove weak convergence, it suffices to show the convergence of the char- 
acteristic function ( |64|) jl|. Indeed, one has 

f hYaMrt) (K,t) = E p (-^(f x (hK) - 1)) h ' r ^ Ep{-t^\K\ a ), (78) 

which completes the proof. □ 

Remark IV. 3. Condition ( 1771) is not void. It is satisfied by all the distributions belonging 
to the basin of attaction of symmetric a-stable laws. Let f a ,x( x ) denote the probability 
density function of a symmetric a-stable law whose characteristic function is 

/^W=e- |s|a , (79) 

then one can immediately see that (1771) holds true. As above, let L a (t) denote the 
symmetric a-stable Levy process. Then, equation (1731) is the characteristic function of 
L a ,p(t) = L a (Np(t)), that is of the symmetric a-stable Levy process subordinated to the 
fractional Poisson process. This remark leads to the conjecture that L a< p(t) is the functional 
limit of Y at p(t), the a-stable compound fractional Poisson process defined by equation (|6T|) 
with the law of the jumps X belonging to the basin of attraction of or coinciding with an 



a-stable law. This conjecture can be found in a paper by Magdziarz and Weron 34| and is 



proved in Meerschaert et al 



and Scheffier (351 ] . 



271 ] using the methods discussed in the book by Meerschaert 
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